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■ Abstract 



The infinite-?/ single impurity Anderson model for rare earth alloys is ex- 



> 

o 

£3 amined with a new set of self-consistent coupled integral equations, which 

(N 

| can be embedded in the large N expansion scheme (N is the local spin de- 

generacy). The finite temperature impurity density of states (DOS) and the 
spin- fluctuation spectra are calculated exactly up to the order 0(1/N ). The 
presented conserving approximation goes well beyond the 1 /TV-approximation 
q [ (NCA) and maintains local Fermi-liquid properties down to very low temper- 



atures. The position of the low lying Abrikosov-Suhl resonance (ASR) in the 



impurity DOS is in accordance with Friedel's sum rule. For N = 2 its shift 
toward the chemical potential, compared to the NCA, can be traced back to 
the influence of the vertex corrections. The width and height of the ASR 
is governed by the universal low temperature energy scale Tk- Temperature 
and degeneracy iV-dependence of the static magnetic susceptibility is found 
in excellent agreement with the Bethe-Ansatz results. Threshold exponents 
of the local propagators are discussed. Resonant level regime (N = 1) and 
intermediate valence regime (|e/| < A) of the model are thoroughly investi- 
gated as a critical test of the quality of the approximation. Some applications 
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to the Anderson lattice model are pointed out. 
PACS numbers: 71.28, 75.20H 
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I. INTRODUCTION 



Recently a major improvement has been reported on the perturbational approach to the 
single impurity Anderson-model (SI AM) with iV-fold degenerate local state in the limit of 
infinitely large Coulomb repulsion iM. The aim of this paper is to present in detail the 
results obtained so far with this so-called Post-NCA theory (PNCA^). This approximation 
contains a resumation of infinite numbers of skeleton diagrams and is exact up to order 
0(1/N 2 ). The skeleton dia grams included physically describe complicated multi-excitation 
processes. 

The paper is organized as follows: in section 2, the large- N expansion is used to derive the 
self-consistent coupled integral equation for the self-energies and the vertex-function. The 
physical processes included will be explained. In section 3 numerically obtained solutions 
for the local propagators will be presented and their threshold exponents related to results 
furnished by other methods. A detailed analysis of the one-particle spectra will be the 
topic of section 4 including a discussion of the Fermi-liquid properties and the impact of the 
vertex-corrections. In section 5, we will develop a theory of the magnetic vertex correction 
exact up to 0(1/N 3 ) and compare the calculated static magnetic susceptibility with the 
Bethe-Ansatz results. Spin-fluctuation spectra will also be covered. The two critical cases, 
the intermediate- valance and the resonant level case, will be thoroughly examined in section 
6. In section 7, a summary and an outlook on further applications will be given. 

II. THEORY 

In direct perturbation theory a well established set of diagrammatical rules has been 
established to evaluate self-energy contributions to the local propagator Pjv/(z)llli. In Fig. 
|l] the lowest order contributions up to V 8 to the self-energy S (z) of propagator Pq(z) of 
the un-occupied state are shown. Contributions (b) and (c) are one-particle irreducible but 
contain self-energy insertions and are already included in the skeleton diagram (a). In order 



3 



V 6 there exist only one skeleton diagram, whereas there are exactly two in order V s . To 
include these higher order corrections, a well-established concept of vertex functions is used. 
In our case, a vertex function is defined by cutting the electron and the local propagator 
line at the uppermost vertex in each skeleton diagram and by summing up all contributions 
of the remaining lower parts of the skeletons. The concept of vertex-function has already 
been successfully applied in an iVCA-theory (Non-Crossing-Approximation) for the finite- 
U Anderson modelS The resulting dimension-less factor A^(x, y) renormalizes the bare 
hybridisation vertex V of an absorbed band electron in an energy dependent way; x is 
the ingoing and y the outgoing energy of the local propagator. The corresponding vertex- 
function A^(x,y) for an emitted band electron is obtained from the vertex function of an 
absorbed electron A^(x, y) by symmetry: A^(x, y) = A^(y, x). Therefore we only calculate 
A m (x, y) = A^(x, y). In a formal way all skeleton diagrams are summed up by using the 
exact vertex-function in: 

£„(*) = V 2 ]T f{e k )A m (z, z + e k )P m {z + e k ) (1) 

km 

S m (z) = V 2 ^2f(-e k )A m (z - e k ,z)P m (z - e k ) . (2) 

k 

Approximations are made by the selection of contributions to A m . The NCA, for example, 
is recovered by setting A m (x, y) = 1. 

In order to derive the integral equation for the vertex function exactly up to 0(1/N 2 ) 
the lower parts of the skeleton diagrams (d) and (e) have to be included in addition to 
skeleton diagram (a). Skeleton (f) can be neglected, since it is of the order 0(1/N 3 ). The 
corresponding generating functional is shown diagrammatically in Fig. [| An infinite number 
of skeletons are resummed by replacing the bare hybridisation vertex by the renormalized 
one V ■ A m (x, y) in each diagram in order to include as many higher order terms as possible. 
This yields a self-consistent equation for the vertex function displayed in Fig. |P. In the case 
of zero magnetic field all P m are equal, and the vertex function is reduced to 
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A(x,y) = A m (x,y) = 

/oo poo 
/ dudvp(u)p(v)f(u)f(—v)G(x,x + u,x + u — v) (3) 
-oo J —oo 

■A(x + u — v,y + u — v)H(y, y — v,y + u — v) . 
Here we have chosen p(e) to be the density-of-states (DOS) of the conduction electrons and 
the auxiliary functions G, H and K 

G(x,x + u,x + u — v) = A(x, x + u)P m (x + u) (4) 

A(x + u — v, x + u)Po(x + u — v) 

H{y } y-u,y + u-v) = P m (y + u - v ) (5) 

■{A{y -v,y + u- v)P (y - v)A(y - v,y) 
-K(y + u-v,y,y + u)) 



K(y + u-v,y,y + u) = N\V\ 2 / dl p(l) f (-1) A(y + u - v - l,y + u - v) (6) 

J — oo 

■Po(y + u — v — l)A(y + u — v — l,y + u — I) 
■P m (y + u- l)A(y -l,y + u-l) 
■P (y-l)A(y-l,y) . 

The negative sign in Eqn.(^) takes into account the odd number of crossing electron lines 
in the skeleton (d) while in diagram (e) an even number of crossing lines can be found. The 
approximation defined by Eqn. (JJ) to (|6|) will be called PNCA^ in the following (oo stands 
for U = oo). 

In the one-particle Green's function F m {iu n ) the vertex correction also comes into play 
natural wayi: 



m a 



1 f dz 

F m (iu n ) = — <p —e'^A^z, z + iu n )P (z)P m (z + iuj n ) . (7) 
Zf Jc 2m 

This equation can e. g. be obtained by cutting one band electron line in the functional $[Pju] 
of Fig. [| Thereby, no ambiguity is left in the analytic expression. The contour C encircles 
the all singularities of the kernel in a counterclockwise fashion. Analytic continuation of the 
Green's function gives 



F m (uj + id) = I dw'P m (u + w') (KA x (w, oo + w')£q(w') 

+^e[P (w')}U x (co,oj + w')) 
dw'P£(w' - uj) {UA y {w' - u, w)f m (tu') (8) 
+ne[P m {w , )}£ Ay {w'-u,w')) . 
We have introduced the real auxiliary functions 

Zm(u) = l —e-^mP M (u + iS) (9) 

itZf 

and the complex auxiliary functions 

%tA x (u), uj + iui n ) = |[A(a> + iS, uj + iui n ) + A{uj — iS, uj + iu) n )} 

QAj.(w, u; + iui n ) = ^t[A(u; + iS, uj + iui n ) — A(uj — iS, uj + iui n )} 

$t,A y (u) — iui n , uj) = ^[A(uj — iu n ,uj + i5) + A{uj — iui ni ui — i5)} 

$sA y (u) — iu n , uj) = ^-[A(c(j — io; n , u; + iS) — A(uj — iu n , uj — i8)] 

UA x ,y) = ~^ e ~ px ^ A x{x,y) 

tA y (x,y) = -^je-^A y (x,y) . 

The functions £m(<^0 introduced above have a simple physical meaning as defect propagator. 
integration over the whole frequency range yields the occupation probability < Xmm > of 
the local state M. The complex functions dt(^s)A x ^ and £,A xy {x,y) enable the evaluation 
of Eqn. (^|) for real frequencies and thus serve purely numerical purposes. 

III. LOCAL PROPAGATORS, THRESHOLD EXPONENTS AND SPECIFIC 

HEAT 

On a cluster of work-stations, we brought a numerical iteration procedure for the above 
system of integral equations to full convergence using dynamically defined logarithmic 
meshes for the threefold integration in (|3]). Besides the modulus \\Pm — P|/ +1 ' ) || (i labels 
the step of iteration), which reaches a value of typical 10~ 16 at the end, the sum rules 

/ ^Pm(z) = 1 ; f ^ m {z) = | V\ 2 £ (< MWJM > [1 - f(e km )\ (11) 

C TV C Til j tm 

<M\{l-nl)\M> f(e km ) 



have been checked to estimate the quality of the numerical calculations. The deviation 
between left and right side is in the range of typically 1-4 % and scales with the inverse 
number of mesh points. For all numerical studies all energies are measured in units of the 
Anderson width A = 7rV 2 Mr {V 2 is the square of the hybridisation matrix element and Mr 
the band DOS at the chemical potential). The featureless symmetric conduction band DOS 
has been chosen to be p c {uj) = 2 r(i 25)iy exp(— (cu/iy) 4 ) to reduce band edge effects. The 
half band width is set to W = 10A. In most cases the temperature will be measured in 
units of the corresponding Kondo-energy 

— • (i2) 

For the investigation of the Kondo-regime e/ = — 3 A and N = 2 is chosen, since the 
largest impact of the vertex function is to be expected for small N. In Fig. |] the spectral 
density poioS) = — SmPo(^ — iS) is displayed in the vicinity of the threshold for three 

7T 

different temperatures T = 0.5, 1, 2T#-. Note the fact that the energy scale has been shifted 
by the threshold energy E s . The transformation Zf = § -^e~ l3z Pm(z) = e~ l3Ea Zf defining a 
renormalized local partition function Zf — § ^ e ~ l3z PM (z) and a local propagator Pm(z) = 
Pm( z + E s ), ensures that the numerically calculated Zf stays of the order 0(1) during the 
variation of temperature. Also the defect propagators £m can be calculated very accurately 
down to very low temperatures. 

With decreasing temperature the threshold behaviour of the propagator P Q in PNCA^ 
exhibits a weaker increase than the corresponding NCA propagator. We regard this as a 
strong hint toward a reduced threshold exponent «o compared with the NCA due to the 
influence of the vertex corrections. On the other hand, the inset indicates quite clearly an 
excellent agreement on the high energy part of the spectrum in both approximations as to be 
expected. The spectral density p m (u), shown in Fig. |]b for the same set of model parameters 
and temperatures, reveals an additional difference between the two approximations. pP NCA 
and P^ NCA develop their common threshold energy E s already at finite temperature while 
in NCA an identical E s for both propagators is found only at T = 0. This gives rise to a 



temperature dependent position of the ASR in the one-particle spectra in NCA, as we will 
see in the next section. 

To obtain a first estimate for threshold exponents we fit the ionic spectrum in the low 
frequency range < u < Tk to a trial function 1im{&) — a (T) ■ u~ aM ^ at different 
temperatures. This procedure has been checked for the NCA and provides the exactly known 
NCA exponents «o = jjzi and a m = within an accuracy of 5%. While no temperature 
dependence is found for the exponent «o = 0.44 ± 0.02, Fig. ||a, the different values of 
a m (T) have been used to to extrapolate a m (0) = 0.27 ± 0.01 for the present parameters 
N = 2,€f = -3 A giving n f = 0.87, Fig. |b. 

There is still no agreement in the literature about the exact exponents for N > 2. The 
proposal of Menge and Muller-Hartmannlll 

«o = ^ a m = ^(2-n f ) (13) 

indicates a dependence of the exponents not only of the degeneracy iV but also on the 
occupation number rif in all regimes. This proposal has been recently backed by a analysis 
of pseudo-boson and fermion propagators provided by the numerical renormalisation groupi. 
On the other side, an analysis of parquet-equations for the Kondo limit of the model (n/ = 1) 
claims that the exponents 

AT 1 _ J_ 

ao = _H = n* (14) 

should be exact in order 0(1/N 2 We did not investigate the threshold behaviour on the 
full scale of the model parameters. Therefore we cannot rule out any of the proposals, even 
though the PNCAoo exponents in the Kondo-regime come very close to Gruneberg's and 
Keiter's results. 

Since we have convinced ourselves that the high energy parts of the spectra remain 
unchanged, we can still use the NCA to obtain specific heat data for a wide temperature 
range. On the other hand, the low energy corrections turn out to be essential for the 
calculation of correlation functions. In Figs. || and |7| origin and scaling behaviour of the 
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two different contributions to the specific heat are clearly demonstrated. In (a) the effective 
hybridisation is tuned leading to a renormalized Anderson width A* and in (b) the position 
of the bare /-electron level has been varied for N = 2. While the first maximum in each 
curve of the specific heat turns out to be an universal function of T/Tk, the position of the 
second rather determined by the /-level energy and its width by A*. Therefore, the first peak 
is generated by spin excitations on the Kondo-scale T ~ Tk, whereas the second contains 
charge excitations from the broadened /-level into unoccupied band states. With increasing 
degeneracy the maximum due to spin excitations is enhanced essentially linearly with 
since their contribution is proportional to the number of channels, as it is already know from 
Bethe-Ansatz calculations^. The charge excitations peak is rather decreased and smeared 
out by broadening of the /-level proportional to iVA. Its position and broadening is in good 
agreement with Bethe-Ansatz0 and numerical renormalisation group calculations^, but the 
height comes out somewhat smaller. This could be due to absence of the doubly occupied 
state for U = oo. Of course, the appearance of the charge excitations is more of academic 
interest for the low temperature behaviour, but may be important for spectroscopic studies. 

IV. ONE-PARTICLE SPECTRA 

The generalized Friedel's sum rule 

nf(T) = - E / fa <W + 7T- E f dzf(z)F m (z)d z E m (z) (15) 

7T m J-oo 2m — Jc 

s » ' 

=k(T) 

relates the local occupation number n/ to the generalized phase shifts 5 m (u) = 
— Smln [— F^{i5)}. k essentially measures the asymmetry of the /-electron spectrum and 
can be interpreted as the negative change of the number of band electrons in the presence of 
the impuritjM In the limit T = 0, Eqn.(15) is used in combination with local Fermi-liquid 



relations to derive the density of states rule: 



^(0,T = 0) = ^-n'(*=tpf>) . (16) 
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This formula predicts a scattering resonance with a maximum height of near the chemical 
potential ((jif — K) ~ 1 ) for N = 2 and T = 0, the ASR, which moves away form fi increasing 
degeneracies N. Fig. [8] shows the one-particle spectra of PNCA^ and NCA for five values 
of temperatures. The most significant differences between both approximations concern 
the position and height of the ASR. The position in NCA is found strongly temperature 
dependent and the height exceeds the unitarity limit of 1/(ttA) already at temperatures 
slightly below Tk- We can trace back the temperature dependence to the mismatch of the 
threshold energies of Pq and P m , which merge only at T = 0. The violation of Eqn. (|16|) 
clearly indicates the importance of the vertex corrections. On the other hand, the PNCA^- 
ASR grows with decreasing temperature at a stable position close to the chemical potential. 
It violates only sightly the density of states sum rule and its agreement with Friedel's sum 
rule ( pT5[) is found to be within 7% as indicated in table |. 

By focusing our attention on Fig. || we gain insight on the impact of vertex corrections on 
the the ASR. The PNCA^ result (solid curve) is compared to the NCA result (dashed curve) 
for the fixed temperature T = 0.5Tk and additionally to a curve where the NCA propagators 
Pm{z) have been used to evaluate the renormalized hybridisation V ■ A(x, y) via Eqn.(^j). In 
the later convergence is achieved after two steps of iteration. The resulting ASR, calculated 
after each step, has been shifted towards the chemical potential, but stays enhanced com- 
pared to the PNCA^-ASR,. The reduction of height of the ASR is clearly connected to the 
modified threshold exponents as is also revealed by an inspection of the iteration procedure. 
This stressed the importance of vertex corrections in the local propagators Pm{z). 



The plots in figure 10 attract attention to the local Fermi-liquid properties as seen in the 



imaginary part of the self-energy of the Green's function: 

3fmE /m (w-i<J) = A + C u uj 2 + C r T 2 . (17) 

From the symmetric Anderson model it is know that C w = A/T^ and C-r = In the 

strongly asymmetric case (U = oo) under consideration here, only the scaling = A/Tj| 
can be found. The quadratic expansion coefficient for the temperature exhibits a rather 
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strong dependence on the degeneracy N as been shown in Fig. [H], originating from the shift 
of the ASR away from the chemical potential, while in the completely symmetric case the 
phase shift 5 m remains 7r/2 independently of N. On the base of our PTVCAoo-results we 
suggest a new analytic investigation of Fermi liquid relations for the asymmetric case. 



V. MAGNETIC SUSCEPTIBILITY AND SPIN-FLUCTUATIONS 

Magnetic excitations are measured by the magnetic susceptibility: 

dz 



-M —M 

z — H z + iv n — H 



X(-Wn) , (18) 



where iv n = ^j 1 is a bosonic Matsubara-frequency. Focusing our attention on the operator 
the local magnetization in the SIAM, M z = gfiBJ2m m '^m,m, we rewrite Xj{iv n ) in analogy 
to F m : 

JSfjj 2 . I r dz 

Xf(iv n ) = tt--^ f —e~ Pz P m {z)P m (z + iu n )r(z, z + iv n ) , (19) 

6 Z Jc Ziri 

defining = j(j + \){g^B) 2 - The new magnetic vertex function V(x,y) formally includes 
all higher order contributions arising from diagrams with crossing band electron lines; in 
NCA T(x,y) = 1. It turns out to be symmetric in its complex energy arguments in order to 
maintain the symmetry xfan) = x{~ iv n ). Using the standard set of diagrammatical rulesi, 
we obtain T(x,y) exactly up to order 0(1/N 3 ) shown in Fig. |ll[ 0(1/N 2 ) contributions 
compensate each other. On the other hand, diagram (c) and (d) in Fig. [TT] transform into 
each other by exchanging the energy arguments. Therefore, we only have to insert 



oo poo 



T c (x,y) = -|U| 4 / / dudvp(u)p(v)f(-l)f(-u)f(v) 

J —oo J —oo 

■Pq(x — u) A(x — u, x) 

■P m (x — u + v) A(x — u, x — u + v) 

■T(x — u + v,y — u + v) 

■P m (y-u + v)j(y + u-v,y,y + u) (20) 
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into the magnetic vertex function of the PA^C/loo-theory 

T(x,y) = l + T c (x,y) + T c (y,x) , (21) 

and to solve Eqn. (p0|) and ( pip self-consistently. Eqn.(^) then leads to a local dynamical 
susceptibility Xfi which is exact in 0{1/N 2 ). 

The quality of the approximation has been checked by comparison of the static magnetic 
susceptibility = lim^o x{ v ) to the exact Bethe-Ansatz results^, Fig. 12. The PNCAoo 



susceptibility shows a surpriselingly good agreement with the exact results, even though 
the Bethe-Ansatz susceptibility has been obtained for the Coqblin-Schrieffer model. In 
particular, the characteristic maximum for degeneracies A^ > 3 has been reproduced in 
PNCAoq. On the other hand in NCA not even a saturation of Xf occurs for 0.1TV < T < Tk 
and N = 2. The deficiencies of the NCA do not show up so strongly at e.g. A^ = 6 or at 
higher temperatures T > Tk, where vertex corrections become less important. 

The scattering function S^) for neutron scattering experiments is linked to the magnetic 
susceptibility by the dissipation-fluctuation theorem: 

s m~t^FM^m • (22) 

The right hand side approaches the spin-fluctuation spectrum a{oj) = x{u)/& for small ui. 



In Fig. [L3| the spin-fluctuation spectra of the PNCA^ and the NCA are compared for three 
different degeneracies AT at a fixed temperature T = 0.27V. The pronounced maxima of 
for N > 3 resemble the maxima in the static susceptibility discussed before. They 
appear on the same energy scale uj(T) m 0.5Tr-. It is also interesting to note that even though 
no maximum is found for N = 2 a. rather broad inelastic peak can be reported in a PNCA (u) 
which is not seen in the NCA spectrum. Of course, this features disappears well above Tk 
for both approximations leaving a single elastic peak at u = in the spectra. The low 
frequency behaviour of Smx(w) can be interpolated by ui /{Yl ieut + oj 2 ). T neut is the neutron 
scattering linewidth which is experimentally defined by the position of the maximum. Here, 
it is convenient to use 
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1 Smx(uj) 



_ lim r_^lZ. (23) 



In Fig. [TJ], the reduced linewidth Y neut /TK is displayed versus temperature for iV = 2, 4, 6 in 



PNCAoo. The linewidth is nearly temperature independent in the local Fermi-liquid regime 
at T < Tr. It develops a weak minimum at T m 0.5T^ which is slightly enhanced by 
increasing degeneracy and behaves like vT for high temperatures. 

VI. CRITICAL EXAMINATION OF THE PNCA^ 

In the last three sections we discussed PNCA^ results obtained for the Kondo regime 
of the model. In this section we will focus our attention on the case N — 1 and the 
mixed- valance regime. For N = 1 all diagrams are of the same order in the sense of an 
1 /iV-expansion scheme. Nevertheless, this classification is questionable here, since there is 
no Kondo effect left although the diagram topology remains unchanged. The one-particle 
Green's function is exactly known: 

F(z) = _^ ; (24) 

Z t f #k z _ ek 

(#k denotes the number of band electron states). In figure |T5| the one-particle spectra of 



the NCA and the PNCA^ are compared with the exact result for two values of 6f. While 
the Pm{ z ) have been determined very accurately, oscillations of the defect propagators 
£(u;) = exp(— f3u)^-QmPM(u) during the iteration procedure cause oscillations in the 
.PiVCMoo-spectra around the exact solution, i.e. convergence is not fully obtained in this 
case. These oscillations contribute about 3% to the residual deviations of the PNCA^. In 
the most critical case, e/ = —A, the TVCA-pathology is strongly pathological near /x, but 
the PNCAoo is much less so. For tf = — 3 A only a very weak pathology at fi — remains, 
while the NCA still produces an unphysical resonance at the chemical potential. 

A rather strange two peak structure is found in the spectrum for the case of the in- 
termediate valence regime 6f = —A and N = 2, shown in Fig. |IT)|a: the first peak nearly 
at /j, = exceeds clearly the density of states rule ( |16|) p(0) = 0.21 for ri/ ps 0.6, which 
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would be the height of the minimum in between both peaks. We suppose that the true 
spectrum will monotonically decrease from the minimum position onwards for decreasing 



energy. The first peak clearly reflects a PiVCAoo-pathology already seen in Fig. |15|a. Despite 
this spurious structure, the overall violation of the the DOS-sum rule is reduced from almost 
60% in NCA down to 15% in PNCA^. In part (b) of the figure, 6f is chosen to be +A. 
Again, only a weak pathology is found in this case, and the spectrum can be reasonably 
fitted using a Lorentzian with a width of 0.95A. This resembles the fact that there is no 
blocking effect in this regime, whereas the local correlations still reduce the spectral weight 
to < X 00 > + < X mm >= 0.87 < 1. 

VII. CONCLUSION AND OUTLOOK 

We have demonstrated that our new FiVCA^-approximation to the Anderson-model 
improves the low energy properties in the Kondo-regime quite remarkably, whereas the sat- 
isfactory behaviour of the old TVCA-approximation at higher energies is maintained. Vertex 
corrections do have in fact a large impact on the low energy excitations: position and height 
of the ASR for N = 2 is shifted towards the chemical potential as consistent with Friedel's 
sum rule. Also the magnetic susceptibility agrees surprisingly well with the Bethe-Ansatz 
results. We plan to apply this method also to the finite-?/ Anderson-model and the ex- 
tended Anderson model&i, which includes a direct exchange interaction between /- and 
conduction electrons. The self-consistent set of equations have been derived alreadyi. Our 
perturbational approach to the Anderson-model has opened the prospect of studying the 
low temperature properties of Heavy Fermion systems in much more detail before. Quite 
generally, the so-called LNCA scheme! can be applied, which provides a very successful 
perturbational treatment of the Anderson lattice. Within this theory based on a picture 
of independent effective sides plus quasi-particle interactions collective effects like super- 
conducting or magnetic ground states can be calculated in a systematic way. In order to 
underline this perspective we present the temperature dependent pseudo-gap formation in 
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the /-spectrum of the Anderson lattice with N = 2 and ej = —3 in Fig. [17]. The temper- 
atures are given in units of the corresponding impurity Kondo temperature which deviates 
moderately from the characteristic temperature of the lattice T* . While the LNCA-scheme 
combined with the local NCA is limited to temperatures larger than Tk due to violations of 
the local Fermi-liquid properties, the so-called L-PNCA^, which uses the PNCA^ to solve 
the effective site, reaches much lower temperatures. Since now the ASR of the impurity 
is located only very slightly above the chemical potential in accordance with Friedel's sum 
rule, the pseudo-gap is formed around uj = in the L-PNCAoc- In a forthcoming publication 
we will investigate the temperature dependent quasi-particle band-structure and the impact 
on the transport coefficients and optical conductivity in greater detail. 
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TABLES 



Friedel's Sum Rule and the NCA-Theory: 



TV TV 
1 1 1 K 


rif 


m 


fares 


0.1 


.8682 


.5863 


.2819 


0.2 


.8682 


.6114 


.2567 


0.3 


.8682 


.6379 


.2303 








(a) 




Friedel's Sum Rule and the PNCA^- Theory: 




1 1 J-K 


Uf 


\ -n 5m 

^ ~K 

m 


fares 


n i 


.8445 


.9023 


0^77 
-.uo 1 / 


0.2 


.8525 


.9140 


-.0615 


0.3 


.8565 


.9207 


-.0641 




TABLE I. Parameter: N = 2, e f = 


-3A, W = 10 


(b) 
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FIGURES 

FIG. 1. Diagrammatic contributions to the self-energy So (z) of the unoccupied propagator Pq 
up to V s . 

FIG. 2. Diagrammatical representation of the self-consistency equation of the vertex function 
A m (x,y) mO(l/N 2 ). 

FIG. 3. Generating functional for the self-energies and the Green's functions up to order 
0(l/N 2 ). The PNCAoc contains additional higher order contributions. 

FIG. 4. Comparison of the spectra of Po(uj) (a) and P m (u) (b) in NCA and PNCAoo at 
T = 0.5, 1.0, 2.0Tk- The increase in height at u> = is correlated with a decrease of tempera- 
ture. Parameters: N = 2, e f = -3A, W = 10 A. 

FIG. 5. Threshold behaviour of the spectrum of Pq(z) (a) and an estimation for the threshold 
exponent a m (b). Parameters as before. 

FIG. 6. Specific heat contribution of the impurity vs temperature T/Tk for (a) different val- 
ues of the effective Anderson width A* = 7r|y*| 2 /9 , ej = — 3A , and (b) different values of e/. 
Parameters: N = 2,e f = -3A , A = 0.1W. 

FIG. 7. Specific heat contribution of the impurity vs temperature T/Aq for (a) different val- 
ues of the effective Anderson width A* = 7r|V*| 2 po, £/ = — 3Ao, and (b) different values of e/. 
Parameters: N = 2,e f = -3A , A = 0.1 W. 

FIG. 8. Temperature dependency of the ASR in (a) PNCA^ and (b) in NCA. Parameters: as 
before. 

FIG. 9. Impact of the vertex correction on the one particle spectrum. Parameters: as before. 

FIG. 10. Minimum of the self-energy versus temperature for iV = 2,4, 6. Parameters: 
€f (N = 2) = -3A, e f {N = 4) = -5A, e f (N = 6) = -6.7A, W = 10A. 
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FIG. 11. Diagrammatic representation of the self-consistency condition of the magnetic vertex 
function T(x, y). 

FIG. 12. Static magnetic susceptibility calculated (a) in NCA, (b) in PNCA^ for the single 
impurity Anderson model and (c) for the Coqblin-Schrieffer models using the Bethe Ansatz. The 
NCA results are normalized with respect to the PNCA^ values at T = 0. Parameters: as in Fig. 10 

FIG. 13. 3mx(w)/w vs u/T K for N = 2,4,6, calculated at T = 0.2T K (a)in PNCA^ and (b) 
in NCA. Parameters: as in Fig. 10 

FIG. 14. Neutron scattering linewidth T neut vs. T/Tk calculated in PNCAoo for different N. 
Parameters: as in Fig. 10 

FIG. 15. One-particle spectra in NCA and PNCA^ compared with the exact solution for (a) 
e f = -A, T = 1/15A and (b) e f = -3A,T = 1/60A. Parameter: N = 1. 

FIG. 16. One-particle spectra in NCA und PNCA^ for (a) e f = -A and (b) e f = +A, the 
intermediate valence regime. Parameters: N = 2,T = ^A. 

FIG. 17. /-spectra of the periodic Anderson model in the vicinity of the chemical potential 
(a) in L-PNCAoo and (b) in LNCA for three temperature T = 0.75, 1.0, 2.0T K . Parameters 
N = 2, e/ = -3A, W = 10A. 
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